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ABSTRACT 


Thus  far,  all  performance  calculations  for  the 
X-Wing  vertical  takeoff  concept  have  used  a  rigid, 
nonflapping  blade  analysis.  At  this  stage  in  the 
concept  development,  the  flapping  degree  of  freedom 
should  be  included  for  increased  accuracy.  Thus,  a 
relatively  simple  set  of  equations  is  developed  to 
relate  rotor  blade  aerodynamic  moments  and  blade 
flapping  response  which  provides  rapid  evaluation  of 
the  blade  dynamics  and  the  resulting  changes  in  aero¬ 
dynamic  loading.  The  relationships  between  the 
aerodynamic  hinge  moments  and  the  blade  response  are 
derived  for  the  first  three  harmonics  of  flapping. 
Changes  in  flap  damping  due  to  flapping  response  are 
accounted  for  in  the  equations,  which  are  solved 
using  standard  matrix  algebra.  The  equations  are 
validated  by  comparison  with  a  more  sophisticated 
analytical  technique.  A  case  is  examined  for  the 
critical  conversion  advance  ratio,  and  the  first- 
order  dynamic  response  is  shown  for  a  range  of  blade 
natural  frequencies.  It  is  shown  that  a  blade 
natural  frequency  higher  than  2-per-rev  is  desired. 
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ADMINISTRATIVE  INFORMATION 

The  work  presented  herein  was  conducted  at  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  (DTNSRDC)  for  the  Defense  Advanced  Research 
Projects  Agency  under  ARPA  order  4238,  DTNSRDC  Work  Unit  1-1690-102. 

For  this  evaluation  all  dimensions  were  either  measured  in  or  converted 
directly  to  U.S.  customary  units.  Hence,  U.S.  customary  units  are  the  primary 
units  in  this  report.  Angular  measurement  is  in  radians  except  where 
degrees  are  specified. 


INTRODUCTION 

The  importance  of  blade  dynamics  to  the  analysis  of  helicopter  performance 
and  control  has  long  been  recognized.  Ongoing  development  of  the  circulation 
control  rotor  technology  and  application  to  the  high-speed  X-Wing  aircraft  have 
concentrated  on  the  unique  aerodynamic  characteristics  of  these  rotor  systems. 
Early  X-Wing  designs  had  very  stiff  blades,  and  an  assumption  of  perfectly  rigid 
blades  was  made  to  simplify  analysis.  However,  interaction  between  the  rotor 
aerodynamics  and  blade  dynamics  does  modify  the  control  phase  angles  for  trim  and 
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the  resultant  hub  moment  responses  to  control  inputs,  especially  on  designs  having 
lower  flapping  stiffness.  Both  items,  therefore,  are  highly  important  to  the 
overall  rotor  and  control  system  design  and  to  the  subsequent  analysis  of 
performance  and  control  sensitivity. 

The  blade  dynamics  model  must  be  easily  incorporated  into  the  existing  rigid- 
blade  aerodynamic  analysis.  An  approach,  therefore,  was  taken  to  develop  a 
separate  blade  dynamics  module.  This  module  uses  aerodynamic  hinge  moments  from 
the  main  analysis,  in  the  form  of  Fourier  series  coefficients,  to  calculate  the 
first  three  harmonics  of  blade  flapping  response.  The  blade  dynamics  model  is  that 
of  a  rigid  blade  with  offset  hinge  and  spring  restraint,  thus  modeling  only  the 
first  natural  flapping  mode  of  the  blade.  Since  the  module  includes  the  first- 
order  dynamic  terms  of  mass,  inertia,  and  centrifugal  force  effects,  these  terms 
are  not  required  in  the  main  program.  The  module  also  includes  aerodynamic  flap 
damping  terms  to  reflect  changes  (up  to  three  harmonics)  that  arise  from  the 
difference  between  the  reference  flapping  and  the  calculated  resultant  flapping 
response.  Aerodynamic  damping  due  to  the  reference  flapping  is  assumed  to  be 
included  in  the  aerodynamic  hinge  moments  provided.  Output  from  the  module  is 
simply  the  Fourier  series  coefficients  of  the  resultant  blade  dynamic  flapping 
response.  These  coefficients  are  then  used  in  the  existing  aerodynamic  analysis 
to  account  for  the  influence  of  blade  flapping  position  and  velocity  on  the  local 
angle  of  attack  and  to  properly  resolve  hub  moments  back  through  the  flapping 
hinge. 

As  a  part  of  this  program  development  effort,  the  dynamics  module  was  used  in 
a  preliminary  study  of  the  potential  influence  of  blade  dynamics  on  the  X-Wing 
rotor.  The  results  of  this  study  are  presented  in  this  report. 

DERIVATION  OF  EQUATIONS 

A  simple  relationship  is  developed  between  blade  hinge  moments  and  blade 
flapping  response,  based  on  the  first  three  harmonics  of  airload  versus  azimuth. 

The  relationship  is  intended  for  use  in  conjunction  with  a  rotor  aerodynamic 
analysis  to  account  for  the  influence  of  blade  dynamics  on  the  rotor  control 
requirements  and  the  proper  control  phase  angles  for  trim  and  maneuver. 
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BASIC  FLAPPING  EQUATION 

The  fundamental  blade  dynamic  relationship  is  given  by: 

IB  +  cB  +  kg  =  M 

e 

where  3  =  flap  angle,  radians 

8,8  =  time  derivatives  of  8 

I  =  mass  moment  of  inertia  of  blade  about  hinge 
c  =  damping  constant 
k  =  hinge  spring  constant 
Mg  =  all  external  moments 

Equation  (1)  describes  the  system  shown  in  Figure  1.  The  three  sources  of 
external  moments  are:  (1)  distributed  lift  forces,  (2)  distributed  centrifugal 
forces,  and  (3)  distributed  blade  weight. 


(1) 


d  l 


Figure  1  -  Offset  Hinge  Dynamics  Model 


EXTERNAL  MOMENTS 

The  external  moments  about  the  hinge  are  expressed  in  integral  form  as: 
R  R 

Mp  *  /  (r-rQ) (d£ /dr)dr  -  f  (r-ro>(3+B  )  rl22  (dm/dr )dr 

i  r  °  P 

o  o 

R 

-  g  /  (r-rQ) (dm/ dr )dr 
r 

o 
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where 


=  hinge  offset  from  centerline  of  rotation,  ft 
n  =  rotational  speed,  rad/sec 
3p  =  precone  angle,  radians 
3  =  flapping  degree  of  freedom,  radians 

Because  d£  and  3  are  both  azimuth  dependent,  the  value  of  M£  is  also  azimuth 
dependent.  The  first  integral  (lift)  is  assumed  to  be  available  from  an 
aerodynamic  analysis  as  a  function  of  azimuth,  M^(<{i).  The  second  integral 
(3  dependent)  and  third  integral  may  be  expressed  in  more  convenient  terms  as: 

Me  =  M^)  -  ft2 (3+3^)  (I  +  room)  -  mgo  (2) 

where  (ip )  =  aerodynamic  hinge  moment  from  distributed  lift  forces 

a  =  center-of-gravity  location  of  hinged  blade  measured  from 
hinge,  ft 

m  =  mass  of  hinged  blade,  slugs 

Combining  Equations  (1)  and  (2)  and  dividing  through  by  the  term  Ift2  yields 
the  following  equation  to  be  solved: 

h  +  tJ1  +  p2s  ‘  ^  -  6f.  (i  *  v"/1'  -  «> 

where  P  =  co  /fi,  first  flap  frequency  ratio 

n 

r  am  , 

<wn/fi)2  =  1  +  -f~  +  W 

y  =  structural  damping  factor 

(=  2t,  two  times  viscous  damping  coefficient) 


FLAP  DAMPING  TERMS 

Local  blade  angle  of  attack  is  changed  by  blade  flapping  response,  as  depicted 
in  Figure  2.  The  amount  of  flap  damping  due  to  the  reference  flapping  condition  at 
the  time  M^(i(/)  was  calculated  is  included  in  the 

new  flapping  solution  will  have  a  different  amount  of  flap  damping,  and  the 
difference  may  be  calculated  from: 


M^(^)  term  of  Equation  (2).  The 


Aa 


tan  1  (w/u) 
-(r-rQ)  AB/u 


(4) 


4 


where  Af3  is  the  difference  between  the  new  solution  flap  rate  and  the  reference 
flapping  condition,  and  u  is  the  local  velocity  vector.  The  aerodynamic  angle  of 
attack  (hence,  load)  and  blade  flapping  are  interdependent. 


Figure  2  -  Change  in  Angle-of-Attack  Due  to  Flapping  Velocity 


This  projected  change  in  flap  damping  hinge  moment  due  to  the  change  in  angle 
of  attack  is  represented  by: 

,R 

=  J  (r-ro)  aAa  q  c  dr  (5) 

o 

Equation  (5)  can  be  evaluated  by  substituting  Equation  (4)  for  Au.  The  u  term 
in  the  denominator  of  Equation  (4)  partially  cancels  the  dynamic  pressure  q  of 
Equation  (5).  The  c  term  in  Equation  (5)  can  be  approximated  by  the  mean 
aerodynamic  chord,  and  the  remaining  u  term  can  be  approximated  by  the  local 
in-plane  velocity.  To  properly  reflect  damping  in  the  reverse  flow  region,  u 
should  be  taken  as  an  absolute  value.  But  to  keep  the  algebra  solvable,  the 
absolute  value  has  been  dropped,  yielding  an  inaccurate  estimate  of  the  damping 
change  in  the  reverse  flow  region.  This  approximation  does  not  influence  the 
aerodynamic  damping  included  in  the  (i^ )  term  of  Equation  (3),  which  correctly 
reflects  all  aerodynamic  terms  from  a  separate  analysis;  therefore,  the  approxi¬ 
mation  does  not  affect  the  converged  flapping  solution.  However,  the  rate  of 
flapping  convergence  for  very  high  advance  ratios  may  be  influenced.  Even  at  an 
advance  ratio  of  one,  the  advancing  blade  damping  is  still  an  order  of  magnitude 
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greater  than  the  retreating  blade  damping.  Thus  the  harmonic  content  of  the 
damping  term  is  only  weakly  influenced  by  the  reverse  flow  region.  Performing 
the  radial  integration  then  gives  a  relationship  between  AM^  and  the  flapping 
solution. 


=  -  (1/2)  YLm(F1  +  pF2  sin  ij,)(3  -  Sref) 


(6) 


where  y 

L 

a 


Lock  number  (pacR4/l) 
lift  curve  slope 


=  1/4  -  2x  /3  +  x  2/2  -  x  4/12 

o  o  o 

=  1/3  -  xo  +  xo2  -  xq3/3 

=  r  /R 
o 

Equation  (6)  represents  a  change  in  the  external  moment  beyond  the  M,  (\{i) 
value  of  Equation  (3).  The  term  AM^/CIS2),  therefore,  must  be  added  to  the  right- 
hand  side  of  Equation  (3)  to  include  the  projected  flap  damping  terms. 


F, 

Fo 


Q2 


+  Y7 


(7) 


The  change  in  the  external  hinge  moment  is  proportional  to  the  difference 
between  the  reference  flapping  and  the  flapping  solution  (B~Bref).  As  dynamic 
and  aerodynamic  solutions  converge,  the  term  (3-B^^)  approaches  zero,  as  does  the 
projected  flap  damping  in  Equation  (6).  Thus,  the  converged  solution  is  unaffected 
by  the  previously  mentioned  approximations  used  in  deriving  the  equation. 


HARMONIC  SOLUTION 

The  simple  harmonic  solution  of  Equation  (7)  is  obtained  in  four  steps: 
1.  The  ^(<10  term  is  assumed  to  be  a  Fourier  Series: 

^(<10  =  Mq  4-  sin  >l>  +  cos  ip 


+  M„  sin  2ip  +  M.  cos  2\p 
2s  2c 


(8) 


+  M.  sin  3iJ;  +  M  cos  3\p 
3s  3c 
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2.  The  solution  for  6  also  is  assumed  to  be  a  Fourier  series: 


8  =  B  + 

o 

eis 

sin 

*  +  B 

lc 

COS 

+ 

e2s 

sin 

2  ip  + 

*2c 

cos 

2>P 

(9) 

+ 

t;3s 

sin 

3ip  + 

63c 

cos 

3* 

3.  The  6  and  8  series  are  obtained  by  differentiating  the  B  series. 

4.  These  series  are  substituted  into  Equation  (7)  and  coefficients  of 
like  terms  are  equated. 

This  procedure  yields  the  following  matrix  formulation  for  the  harmonic  solution 
of  the  flapping  coefficients  in  Equation  (9)  in  terms  of  the  excitation  >.ffi- 
cients  from  Equation  (8) : 


A{B}  =  {M} 


(10) 


where 


{B}  = 


Is 


'lc 

K> 


"2c 


3s 


l.B3cj 


column  array  of  flapping  solution  harmonics 


and 
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r 


m  r  a  , 

f  M‘o  -  Sp  +  -r->  -  ir  -  {  ylf2p  su 

M‘is-K  <Vtc  +  F2“e2.> 

1  +  F2UB2c> 


*  \ 


lc  2  L  lMls 


{M}  =  {  Ml2s  -  I  <2F1B2c  -  2  V6L  +  I  VV  / 

Ml2c  "  2  YL  (_2F1B2s  ”  2  F2uB1c  +  2  F2VJg3c) 
M1,  ~  i  Yt  (3F  ft*  -  F,u0*  ) 


M 


3s  2  'L  w  lw3c  2  2s 


7  Y,  (-3F.fi  -  F„pg  ) 


V.  3c  2  'L 


1m3s  2hm2cj 


column  array  of  external  moment  harmonics 


Here 


M1  .M1 
ns  nc 


*  * 

6  ,g 
ns  nc 


coefficients  of  the  sin(ni p)  and  cos(ni|;)  terms,  respectively, 
from  Equation  (8)  divided  by  (IQ2) 

coefficients  of  the  harmonic  terms  of  the  reference  flapping 
condition 

harmonic  number  (1,  2,  or  3) 


Also, 


”  all 

0 

al3  ° 

0 

0 

0 

0 

a22 

a23  a24 

0 

0 

0 

0 

a32 

a33  ° 

a35 

0 

0 

0 

a42 

0  a44 

a43 

a46 

0 

0 

0 

a53  a54 

a55 

0 

a57 

0 

0 

0  a64 

0 

a66 

a67 

O 

_ ) 

0 

0  0 

a75 

a76 

a77  -> 

8 


where 


all 

a22 

a33 

a44 

a55 

a66 


P2, 

a13 

=  -YLF2M/4 

P2-1, 

a23 

=  -(YP+ylF1/2), 

a24 

P2-1, 

a32 

=  (yP+Y^/2), 

a35 

P2-4, 

a45 

=  -2(yP+YlF1/2), 

a42 

P2-4, 

a54 

=  +2(YP-PylF1/2), 

a53 

P2-9, 

a67 

=  -3(yP+YlF1/2), 

a64 

P2-9, 

a76 

=  +3(YP+YlF1/2), 

a75 

=  -YLF2p/2 

=  -YLF2u/2 

=  YlF2u/4, 

a46  =  "3ylVM 

=  YLF2P/4, 

a57  =  "W'4 

=  YlF2p/2 

=  YlF2p/2 

Equation  (10)  is  then  solved  by  standard  matrix  methods  to  yield  the  solution 
vector  { B } -  Note  that  the  first  equation  represented  in  the  matrix  set  is  that  of 
the  flap  coning  angle,  This  equation  may  be  decoupled  from  the  rest,  leaving 

only  a  6  by  6  matrix,  for  some  savings  in  computation  time.  The  term  can  then 
be  written  as  a  separate  algebraic  expression  in  terms  of  to  which  it  is 

coupled. 


XWDYN  COMPUTER  PROGRAM 


DESCRIPTION 

The  harmonic  solution  is  programmed  in  Basic  2.0  to  run  on  an  HP-9836  desktop 
computer.  This  program  has  fixed  inputs  to  represent  (1)  the  blade  dynamic 
characteristics,  (2)  the  aerodynamic  hinge  moments  versus  azimuth  as  a  Fourier 
series,  and  (3)  the  initial  blade  flapping  angle  versus  azimuth  as  a  Fourier 
series.  The  solution  is  a  direct  application  of  Equation  (10). 

VALIDATION 

The  X-Wing  dynamics  module  XWDYN  was  validated  against  a  rotor  blade  dynamics 
solution  from  the  Circulation  Control  Rotor  Performance  Program,  CCRPERF.  The 
CCRPERF  program  solves  blade  dynamics  by  azimuthal  integration  of  the  coupled 
flap-lag  degrees  of  freedom  until  a  harmonic  solution  is  obtained,  i.e.,  dynamic 


parameters  at  the  end  of  a  revolution  are  nearly  the  same  as  those  at  the 
beginning  of  the  revolution. 

There  are  several  significant  differences  between  the  two  dynamics  models. 
First,  XWDYN  is  restricted  to  only  the  first  three  harmonics;  whereas,  CCRPERF  is 
not  restricted.  Second,  XWDYN  is  for  flapping  only;  whereas,  CCRPERF  yields  a 
solution  to  the  coupled  two  degree-of-f reedom  flap-lag  problem.  Third,  the 
CCRPERF  integration  assumes  that  aerodynamics  are  constant  over  the  azimuth 
interval;  whereas,  XWDYN  takes  a  continuous  Fourier  series  representation  for 
aerodynamics.  Finally,  the  CCRPERF  model  accounts  for  kinematic  pitch-flap 
coupling  (delta  3),  which  introduces  changes  in  the  effective  stiffness  of  the 
blade.  As  an  isolated  routine,  XWDYN  does  not  account  for  pitch-flap  coupling, 
but  includes  that  influence  as  reflected  in  the  input  aerodynamic  hinge  moments. 

A  validation  case  was  run  on  the  CCRPERF  program.  The  trim  condition  was 
chosen  to  have  non-zero  hub  moments  to  ensure  some  blade  flapping  response. 
Aerodynamic  hinge  moments  and  corresponding  blade  flapping  conditions  were  taken 
from  the  converged  CCRPERF  solution  and  used  as  input  to  the  XWDYN  routine. 
Calculated  flapping  response  from  XWDYN  must  match  that  from  CCRPERF  to  validate 
the  XWDYN  program. 

Typical  output  from  the  computer  runs  is  shown  in  Figure  3.  The  aerodynamic 
hinge  moments  and  reference  flapping  condition  are  specified  input  parameters. 

The  XWDYN  solution  results  in  a  1/rev  flapping  magnitude  within  2.7  percent,  and 
a  phase  angle  within  1.5  deg,  of  the  CCRPERF  1/rev  flapping.  Considering  the 
above-mentioned  differences  between  the  two  dynamic  analyses,  the  validation 
comparison  shows  good  agreement  between  the  CCRPERF  reference  flapping  condition 
and  the  XWDYN  flapping  solution. 

DYNAMIC  RESPONSE  STUDY 

The  XWDYN  module  was  applied  to  an  X-Wing  case  to  evaluate  the  module's 
influence  on  the  "rigid"  blade  aerodynamics.  These  results  are  significant  in 
the  relative  changes  from  case  to  case  as  well  as  in  the  differences  from  the 
reference  "rigid"  case.  However,  accurate  solutions  for  the  amount  of  flapping 
require  that  XWDYN  be  coupled  with  a  complete  aerodynamics  model  to  provide  an 
iterative  solution  of  flapping,  aerodynamic  loads,  and  hinge  moments. 
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Figure  3  -  Inputs  and  Results  for  Validation  Case 


PHYSICAL  CHARACTERISTICS 

The  X-Wing  case  chosen  was  that  of  a  full-scale,  50-f t-diameter  rotor 
operating  at  0.7  advance  ratio.  The  aerodynamic  hinge  moments  were  calculated 
at  DTNSRDC  about  a  34-percent  offset  hinge.  The  rotor  was  intentionally  out-of¬ 
trim  to  provide  a  non-zero  flapping  solution.  Pertinent  parameters  for  the  case 
are: 


Thrust 
Hub  Roll 
Hub  Pitch 
Shaft  Angle 
Collective  Pitch 
Cyclic  Pressure 


12,145  lb 
3,666  ft-lb 
11,253  ft-lb 
+6  deg 
-5  deg 

(in  terms  of  pressure  ratio) 


Sin(ij;)  -0.25 

Cos(i(j)  0.00 
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From  these  parameters,  it  is  evident  that  the  positive  hub  roll  is  due  to  the 
combination  of  negative  sin(i^)  control  input  and  collective  pitch  setting; 
whereas,  the  positive  hub  pitch  is  a  result  of  positive  shaft  angle  and  zero 
cos(ijj)  control  input.  These  observations  for  the  rigid  rotor  case  are  important 
for  a  later  interpretation  of  the  dynamic  responses. 

Dynamic  characteristics  are  those  of  a  typical,  50  ft-diameter  X-Wing  design. 
Mass  and  inertia  terms  were  estimated  from  structural  requirements  based  on  the 
34-percent  hinge  offset.  Typical  output  for  the  design  flapping  frequency  case 
P  =  1.577  is  shown  in  Figure  4.  The  most  obvious  feature  of  the  graph  is  the 
dominant  2/rev  flapping  content.  This  is  consistent  with  the  large  2/rev  content 
from  the  aerodynamic  hinge  moment.  The  curves  in  Figure  4  illustrate  the 
importance  of  including  the  projected  aerodynamic  flap  damping  terms  in  the 
solution.  The  curve  without  projected  flap  damping  shows  a  much  larger  2/rev 
content  than  the  curve  with  projected  flap  damping.  The  projected  flap  damping 
is  attributable  to  the  change  in  flapping  from  the  reference  flapping  state.  In 
the  sample  case,  these  X-Wing  aerodynamic  loads  are  from  a  nondynamic  model  with 
zero  reference  flapping.  Thus,  all  of  the  aerodynamic  flap  damping  present  in  the 
solution  shown  in  Figure  4  was  added  by  the  XWDYN  program. 
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Figure  4  -  Flapping  Response  for  X-Witig  Case 
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The  potential  influence  of  flapping  motion  may  be  estimated  from  the 
corresponding  change  in  blade  element  angle  of  attack.  An  approximation  of  the 
change  in  blade  tip  angle  of  attack  was  made  using  Equation  (4)  and  is  shown 
versus  rotor  azimuth  angle  in  Figure  5.  For  i|i  =  0  to  180  deg,  the  change  is 
moderate,  oscillating  about  +1  deg;  however,  much  larger  changes  in  angle  of 
attack  are  evident  at  other  azimuth  angles.  These  changes  range  from  -4  deg  at 
ip  =  220  to  +9  deg  at  ip  =  290  deg,  which  suggests  that  flapping  motion,  even  though 
small,  strongly  influences  aerodynamic  loads. 


AZIMUTH  ANGLE  (deg) 

Figure  5  -  Approximate  Change  in  Blade  Tip  Angle  of  Attack 
FLAP  FREQUENCY  SENSITIVITY 

The  sample  solution  described  was  performed  over  a  range  of  blade  flap 
frequency  ratios  to  examine  the  degree  of  variation  in  blade  response.  All  other 
values  were  held  constant  during  the  flap  frequency  changes.  Input  aerodynamics 
were  originally  from  a  nondynamic  rotor  analysis  and,  therefore,  should  apply 
equally  well  to  all  cases.  The  blade  hinge  offset  and  the  mass  and  inertia 
properties  were  held  constant  as  originally  calculated  for  the  flap  frequency 
ratio  case;  see  Figure  4.  Thus,  the  frequency  increases  must  be  considered  a 
result  of  increased  hinge  spring  stiffness. 
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Results  of  this  study  show  that  flap  frequency  has  a  profound  influence  on 
the  character  of  flapping  response.  Flap  frequency  also  can  significantly  alter 
the  aerodynamic  load  distributions  from  those  of  a  "rigid"  blade,  nondynamic  model 
(at  least  for  the  u  =  0.7  case  with  non-zero  hub  moments  considered  here).  The 
1/rev  content  of  net  hinge  moment  is  also  modified  considerably.  However,  the  net 
effect  of  flapping  motion  on  hub  trim  moments  must  be  determined  by  the  more 
complete  math  model. 

One  example  of  the  flap  frequency  influence  on  loads  occurs  as  the  blade 
frequency  crosses  the  2/rev  boundary.  A  simple,  single  degree-of-freedom  (mass, 
spring,  damper)  system  will  change  the  phase  lag  of  its  response  when  the 
excitation  frequency  goes  from  below  system  natural  frequency  to  above  system 
natural  frequency.  This  also  holds  true  for  the  rotor  blade  response  as  its  first 
natural  flapping  frequency  is  moved  (analytically)  about.  Figure  6  shows  blade 
flapping  response  for  flapping  frequency  ratios  of  1.8P  and  2.2P.  The  undamped 
flapping  response  (dash  curves)  is  dominantly  a  cos(2ijj)  content,  which  dramatically 
changes  sign  between  the  1.8P  case  and  the  2.2P  case.  The  damped  response,  how¬ 
ever,  behaves  in  a  much  more  benign  manner,  with  no  change  in  2/rev  magnitude  and 
only  a  25-deg  change  in  2/rev  phase  angle.  Here  again  is  vivid  evidence  of  the 
important  role  damping  plays  in  predicting  blade  dynamic  response  during  conversion 
and  hence  its  importance  in  loads  and  vibrations.  The  stiffness  increase  also 
yields  a  substantial  (5:1)  reduction  in  the  1/rev  flapping  magnitude  (not  readily 
visable  from  the  figure).  These  combined  influences  reduce  the  1/rev  net  hinge 
moment  of  the  2.2P  case  to  only  40  percent  of  that  for  the  1.8P  case. 

-  WITH 

PROJECTED 
FLAP 
DAMPING 

- WITHOUT 

PROJECTED 
FLAP 
DAMPING 

0  90  180  270  360  “  ~0  90  180  270  360 

AZIMUTH  ANGLE  (dagl 

Figure  6a  -  P  =  1.8  Figure  6b  -  P  =  2.2 

Figure  6  -  Comparison  of  Flapping  Response  for  P=1.8  and  P=2.2 
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The  variation  in  overall  flapping  response  with  changing  frequency  ratio  is 
summarized  in  Figure  7.  Flap  angle  magnitude  and  azimuth  angle  at  maximum  flap 
are  shown  for  each  of  the  first  three  harmonics  over  a  range  of  frequency  ratios. 
The  curves  are  extended  to  a  frequency  ratio  of  10,  where  the  "rigid"  blade 
solution  is  plotted.  The  dominate  2/rev  response  rapidly  diminishes  as  the  blade 
natural  frequency  ratio  extends  beyond  2.3/rev.  The  limited  amplitude  response 
at  resonance  (2-P  crossing)  reflects  a  significant  (flap)  damping  content.  The 
phase  angle  of  the  2/rev  response  shows  a  gradual  shift,  leveling  out  as  the  blade 
flap  frequency  ratio  approaches  3/rev.  The  3/rev  flapping  magnitude  maintains  a 
moderate  level  until  the  blade  flap  frequency  ratio  goes  beyond  3/rev.  This 
resonant  response  near  3/rev  is  somewhat  academic  because  the  3P  crossing  will 
only  occur  at  a  reduced  rotor  rpm,  and  hence  at  a  higher  advance  ratio  than  that 
of  these  cases.  Thus,  the  aerodynamic  loads  at  p  =  0.7  do  not  apply.  The  case 
is  still  useful,  however,  as  a  check  on  the  validity  of  the  approach. 

The  1  P  flap  resonse  shown  in  Figure  7  is  indicative  of  control  response  in 
that  the  magnitude  of  the  1  P  net  hinge  moments  passed  to  the  hub  are  proportional 
to  the  1  P  flapping  magnitude.  This  results  in  a  steady  moment  applied  to  the 
airframe.  The  1  P  flapping  content  rapidly  diminishes  as  flap  frequency  ratio 
increases  to  a  ratio  of  about  2.2/rev.  This  is  followed  by  a  90-deg  phase  shift 
(at  nearly  constant  amplitude)  as  the  flap  frequency  ratio  increases  to  about 
2.6/ rev. 

NET  HINGE  MOMENTS 

The  net  hinge  moment  at  each  azimuth  is  calculated  from  the  product  of  hinge 
spring  stiffness  and  flap  angle  deflection  for  each  azimuth.  Thus,  the  1/rev 
content  of  the  net  hinge  moment  is  obtained  by  direct  multiplication  of  the  spring 
stiffness  by  the  1/rev  flapping  content.  The  relation  that  this  net  moment  has  to 
total  hub  moment  varies  with  shear  load  at  the  hinge  and  any  aerodynamic  loads 
inboard  of  the  hinge.  Still,  it  is  of  interest  to  observe  the  net  hinge  moment  as 
it  varies  with  flap  frequency  ratio,  because  the  net  hinge  moment  is  a  major 
portion  of  the  total  hub  moment  and  vehicle  trim  moments. 
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Figure  8  shows  the  sine  and  cosine  components  of  1/rev  net  hinge  moment 
versus  flap  frequency  ratio.  Note  that  net  hinge  moment  must  approach  zero  as  the 
spring  stiffness  approaches  zero,  regardless  of  the  flap  angle  magnitude.  The 
"rigid"  blade  values  are  plotted  at  a  frequency  ratio  of  10.  As  stated,  the 
negative  sini|*  pneumatic  control  input  (and  collective  pitch  setting)  creates  a 
negative  sini|;  aerodynamic  hinge  moment.  As  shown  in  Figure  8,  the  sinij/  component 
of  net  hinge  moment  is  always  of  lower  magnitude  than  that  of  the  rigid  blade. 

The  lowest  frequency  ratio  of  1.577  yields  about  a  2:1  reduction  in  net  hinge 
moment,  suggesting  some  reduction  in  roll  control  effectiveness.  For  the  "rigid" 
blade  case,  there  is  no  cosi|/  pneumatic  control  input,  but  there  is  a  negative  cost); 
aerodynamic  hinge  moment  arising  from  the  positive  shaft  angle.  This  yields  a 
strong  pitch-up  hub  moment.  The  costjj  component  of  net  hinge  moment  varies  con¬ 
siderably  with  flap  frequency  ratio. 


Figure  8  -  Variation  of  One/Rev  Net  Hinge  Moments 
with  Flap  Frequency 
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Net  hinge  moment  is  shown  in  Figure  9  as  a  cross  plot  of  the  cosp  versus  slnip 
components  with  flap  frequency  as  a  parameter.  The  sint|/  and  cosij/  net  hinge  moments 
have  a  consistent  trend  toward  the  "rigid"  blade  values  as  frequency  is  increased. 
Figure  9  shows  that  the  most  simp  response  for  a  simp  control  input  and  the  least 
cosip  response  for  no  cosi|/  control  input  occur  in  the  2.6/rev  to  2.9/rev  range  of 
flap  frequency  ratio.  Conversely,  the  minimum  sin^  response  and  the  most  netative 
cost]/  response  tend  to  occur  near  the  2/rev  frequency  ratio.  These  responses  are 
about  -90  deg  out  of  phase  with  the  pneumatic  control  input  signal.  The  best  1/rev 
control  response,  therefore,  is  judged  to  occur  for  flap  frequency  ratios  in  the 
2.6/rev  to  2.9/rev  range. 


NET  HINGE  MOMENT  -  SINE  TERM 

Figure  9  -  Cross  Plot  of  Net  Hinge  Moment  Components 

APPLICATION  TO  ROTOR  PERFORMANCE  ANALYSIS 


HUB  MOMENTS 

To  evaluate  rotor  hub  moments  properly,  the  flapping  solution  must  be  used. 
Flapping  motion  changes  the  azimuth  angle  where  the  combined  aerodynamic  and 
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inertia  loads  are  passed  back  through  the  hinge  to  the  hub.  Referring  to  Figure  1, 
the  hub  moment  is : 


r  R 

HM(iJj)  =  J°  r(d!l/dr)dr  +  r  /  (d£/dr)dr  +  k  B(^) 
o  r 

o 


r  mo  B  ( ijy ) 
o 


(ID 


The  first  term  is  the  hub  moment  due  to  aerodynamic  lift  inboard  of  the  hinge,  if 
any.  The  second  term  is  the  hub  moment  from  the  shear  force  acting  at  the  hinge. 
The  third  term  is  the  moment  passed  across  the  hinge  spring;  the  fourth  term  is 
the  inertial  shear  at  the  hinge.  These  terms  implicitly  account  for  all  of  the 
inertia,  lift,  and  centrifugal  forces  acting  on  the  blade,  as  reflected  in  the 
flapping  angle  solution  B.  If  the  hinge  is  a  pure  pin  (no  spring),  then  k=0,  and 
no  hinge  moment  exists.  Thus,  flapping  changes  the  hub  moment  by  modifying  the 
aerodynamic  load  through  the  8  influence  and  by  introducing  dynamic  terms. 


SOLUTION 

The  XWDYN  dynamics  module  may  be  added  to  a  rotor  aerodynamics  program  in  a 
general  (modular)  way,  thus  allowing  for  more  elaborate  flapping  solutions  in  the 
future.  Figure  10  outlines  a  technique  where  the  specific  airload-to-f lap  relation 
is  accounted  for  in  the  dynamics  module.  As  described,  the  module  basically  would 
be  that  of  Equation  (10).  Other  necessary  modifications  to  the  main  program,  for 
this  or  any  other  dynamics  module,  include: 

•  Providing  for  dynamic  characteristics  in  main  I/O  routine 

•  Calculating  the  aerodynamic  hinge  moment  versus  azimuth,  M^(i|j) 

•  Using  Fourier  series  for 

•  Adding  Aa  term  due  to  8  in  aerodynamic  loads 

•  Adding  a  dynamic  equilibrium  check 

•  Changing  hub  moment  calculations  to  reflect  the  dynamic  module 

The  advantage  in  using  this  solution  is  that  it  is  easily  added  to  an  existing 
main  program.  While  more  elaborate  techniques  may  be  used,  this  technique  has  the 
least  impact  on  the  original  main  program.  The  modular  approach  allows  future 
improvements  to  the  flapping  representation. 


19 


* 

** 

*** 


** 


Existing 

Requires  Modification 
New  Module 


Figure  10  -  Modular  Inclusion  of  Rotor  Blade  Dynamics 


SUMMARY 

A  relatively  simple  relationship  has  been  derived  to  allow  calculation  of  the 
first  three  harmonics  of  blade  flapping  response  to  the  specified  first  three 
harmonics  of  aerodynamic  hinge  moments.  This  relationship  was  developed  as  a 
computer  program  dynamics  module,  XWDYN.  This  module  can  be  easily  incorporated 
into  existing  rotor  aerodynamic  analyses  to  account  for  the  first-order  influence 
of  blade  flapping  response.  The  XWDYN  dynamics  module  has  been  validated  by  compar¬ 
ing  results  with  those  from  another  analytical  technique. 

A  brief  study  was  made  to  evaluate  the  potential  influence  of  blade  dynamics  on 
the  solution  of  an  X-Wing  conversion  case.  Results  show  high  2/rev  blade  response 
for  the  particular  case  examined  (p  =  0.7).  The  magnitude  and  phase  angle  of 
response  change  dramatically  with  changes  in  the  first  flap  frequency  of  the  blade. 
In  particular,  this  response  shows  a  strong  influence  on  both  the  2/rev  flapping 
and  on  the  first  harmonic  of  net  hinge  moment,  suggesting  a  frequency  influence  on 
rotor  cyclic  control  effectiveness.  This  effect  should  be  examined  further  by 
incorporating  the  XWDYN  module  into  a  rotor  aerodynamics  routine  for  a  complete 
iterative  solution. 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  TECH¬ 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDLESS  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS.  A  SEMI  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  A  PRELIM¬ 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  INTEREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN¬ 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMERICAL  CODE  OF  THE  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE-BY-CASE 


